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BEST INDICATOR ADAPTIVE 
FORECASTING METHOD 

DESCRIPTION 
BACKGROUND OF THE INVENTION 

5 Field of the Invention 

The present invention generally relates to a computer implemented 
method of forecasting product demand and, more particularly, to a unifying 
forecasting framework called "Best Indicator Adaptive" or BIA method which 
encompasses many individual forecasting systems, each making use of single, 
10 double or triple indicators while sharing a central common theoretical 

foundation as well as a global framework and methodology uniting all the 
indicators together to produce a final optimum forecast. 

Background Description 

Traditional time series statistical forecasting makes use of only 
15 demand history, that is, demand in time periods in the past, and project to the 

future, assuming that patterns in the past will repeat in the future. Although 
some have made attempts to use orders (load) of current time period in 
making a forecast, none have made use of a variety of information and 
indicators all related to demand in the current time period such as load, ship, 

20 CA (customer accept) history, and exploit the relationships among them in 

/ 

making a forecast. Neither has anyone in the past made use of the aggregated 
pattern in the dates for the orders to be fulfilled in the future to make a better 
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forecast. Finally, none has a process to adaptively choose the best model 
among those just described to come up with the final optimum forecast. 



SUMMARY OF THE INVENTION 



It is therefore an object of the present invention to provide a unifying 
5 forecasting framework which encompasses many individual forecasting 

systems, each making use of single, double and triple indicators. 

According to the invention, the system makes use of four sources of 
information, creating seven different forecasting models. The adaptive 
optimization finally makes use of these seven models to produce a final 
10 forecast. The invention significantly reduces the forecast error for any given 

individual indicator or forecasting subsystem. 

The four sources of information or indicators are the following: 

1 . Load or total order (L)\ 

2. Ship (5); 

15 3 . C A Quarterly history (CA to ); and 

4. CRAD (customer requested date) or RSD (requested ship date for the 
load or orders. 

In the forecasting framework according to the invention, a plurality of 
forecasting subsystems are incorporated, but only one among the plurality 

20 makes use of the information in the past only. In a specific implementation of 

the invention, seven forecasting subsystems are incorporated. All these seven 
forecasting methods share the same central fundamental theoretical 
foundations while each maintains its own uniqueness. A unique capability of 
the invention is the optimization framework making use of all the seven 

25 indicators. This novel and unique capability significantly reduces the forecast 

error for any given individual indicator or forecasting subsystem. 
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BRIEF DESCRIPTION OF THE DRAWINGS 



The foregoing and other objects, aspects and advantages will be better 
understood from the following detailed description of a preferred embodiment 
of the invention with reference to the drawings, in which: 

Figure 1 is a block diagram showing the overall system in which the 
invention is implemented; 

Figure 2 is a system flow diagram illustrating the process implemented 
by the invention; 

Figure 3 is a flow diagram showing the process for the function 
forecast generation from load (X); 

Figure 4 is a flow diagram showing the process for the function 
forecast generation from ship (5); 

Figures 5A and 5B, taken together, are a flow diagram showing the 
process for the function forecast generation from load and ship (LS); 

Figures 6A and 6B, taken together, are a flow diagram showing the 
process for the function CA LCKAD forecast generation; 

Figure 7 is a flow diagram showing the process for the function of 
adaptive optimization; 

Figure 8 is a graph showing the relationship between the ratios 
Load/CA and Ship/Load and a non-BIA forecast; 

Figure 9 is a graph showing the comparison between the BIA's 
forecast making use of the relationship between the ratios Load/CA and 
Ship/Load and a non-BIA forecast; 

Figure 10 is a graph showing the relationship between the forecast 
adjustment and the signal-to-noise ratio (SNR); 

Figure 1 1 is a graph, similar to Figure 10, but showing that the point to 
be forecasted is not included in the fitting of the regression model; and 
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Figure 12 is a block diagram of the system using CRAD for enhancing 
forecast accuracy with Load. 

DETAILED DESCRIPTION OF A PREFERRED 
EMBODIMENT OF THE INVENTION 

5 Referring now to the drawings, and more particularly to Figure 1, there 

is shown the overall system diagram. There are two kinds of data. One is 
historical data 10 and the other, current data 1 1 . For example, for load, 
historical data 10 refers to the load that occurs in the previous time periods or 
quarters. Current data 1 1, on the other hand, refers to the order that is for this 

10 present quarter where the forecast is to be made. BIA 12 takes the data input 

and produces the forecast output 13. Such output is fed to the supply decision 
high level executive meeting 14, called SOP meeting (Sales and Operations 
Meeting). The executives take into consideration the forecast 13 made by BIA 
12 as well as other business information such as profit margin, product 

15 commonality, marketing campaigns, to make the final supply decision 15. 

Such decision drives the factory 16 and the supply chain in producing the 
products (e.g., computers 17) to satisfy the customer demand, which in turns 
brings back revenue 1 8 and profit to the business unit. 

Without a good forecast, the executives would be blind in making 

20 supply decision. Supplying too much information results in scraps and 

excessive inventory. Supplying too little information results in lost revenue 
and customer satisfaction. Either way is detrimental to the vitality of the 
business. Making accurate forecast is crucial to the success of the business. 

Figure 2 shows the system functional diagram for the BIA 12 in Figure 

25 1 . The data are classified into four sources (though the first two, Load and 

CRAD, are regarded as load related information). Load 201, Ship 202 and CA 
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quarterly history 203 are each used to create respective forecasts 204, 205 and 
206 based on the single indicator source information alone. By modeling the 
ratio of quarter-to-date load to quarter CA actual as a random variable with 
gamma distribution, the CA becomes a variable with generalized gamma 
5 distribution whose mean and sigma can be easily computed from the sample 

mean and sigma of the Load-to-CA ratio. The outputs generated by forecasts 
204, 205 and 206 respectively are CA L 207, CA S 208 and CA hisl 209. Load 
201, Ship 202 and CA history 203 are used as a trio to support the generation 
of the forecast 210 called CA^ 211. Because the Load to CA ratio exhibits 

10 significant uncertain and large sigma for its distribution from time to time, the 

estimator for CA^ 21 1 makes use of a unique property not previously known 
anywhere that the ratio of Load-to-CA relates very well to the Ship-to-Load 
ratio. By estimating the functional relationship and the parameters relating 
these two ratios, BIA can predict the load-to-CA ratio with much less 

15 volatility or sigma, making use of the current Ship-to-Load ratio and the 

functional relationship relating the two ratios. Once the load-to-CA ratios are 
estimated with higher certainty, the final CA forecast can be produced with 
higher certainty also. 

The next feature of the invention has to do with how the forecasts 

20 made by the above mentioned methods can be refined using the CRAD 

information. In Figure 2, CA L 207, CA 5 208 and CA^ 21 1 are used for 
respective forecast generation and CRAD 212, 213 and 214. These, in turn, 
are used to refine the previous three forecasts as CA £CRAD 215, CA 5CRAD 216 
and CA^cj^ 217, respectively. The key idea behind it is sometimes the order 

25 could be built up artificially (large quantity but weak). When such a situation 

occurs, generally more orders are requested to be shipped later than usual. By 
taking the log of the ratio of the mean over sigma for the distribution of the 
orders according to the dates the orders are to be shipped, such log signal-to- 
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noise ratio (SNR) becomes a good predictor of the relative strengths of the 
demand. It predicts the relative offset of the forecast made by the previous 
methods (CA^, CA 5 , CA^) to create new forecasts (C A £ CRAD , CA 5CRAD , 
CA^cr^). Such a predictor proves to be very valuable in improving forecast 
5 accuracy when applied to demand forecasting. 

The functional box 218 called "Adaptive Optimization" makes use of 
the seven forecast models created prior to that stage (CA L , CA S , CA^, 
CA L CRAD , CA 5CRAD , CA^cr^, and CA his ) and create a final optimum forecast 
219. There are several keys to this function. One is that it picks the best 

10 forecast model specific to the geography and product group the forecast is to 

be applied to. This is very crucial to the success of the function. The second 
key is that it eliminates candidates depending on known properties according 
to how long the historical data are available and whether the time of forecast is 
early in the quarter, late in the quarter or not. 

15 Forecast generation from Load 204 is performed as shown in the flow 

diagram of Figure 3. In function block 31, the ratios of load to CA are formed 
for each of the historical quarter / at the same week j. Then in function block 
32, the sample mean and sigma of the collection of history of the ratios are 
computed. If the load-to-CA ratio is modeled by a gamma distribution with 

20 parameters a, p, then the distribution for the final forecasted demand becomes 

a generalized gamma distribution. So, in function block 33, the two 
parameters a, p are estimated, and in function block 34, the mean and sigma 
of the generalized gamma distribution are estimated for the final forecasted 
CA demand. 

25 Forecast generation from Ship 205 is performed as shown in the flow 

diagram of Figure 4. This process is very similar to the last function in Figure 
3. The only difference lies in the L fi in the numerator for the ratio in function 
block 31 (now function block 41) and the quarter-to-date current week load Lj 
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used in function block 34. They are changed to quarter-to-date ship S y , and the 
current week ship S y to compute the sample mean and sigma of the collection 
of history of the ratios. 

The forecast generation from Load/Ship (LS) 210 is performed as 
shown in the flow diagram of Figures 5 A and 5B. This function in the 
invention exploits the relationship between the ratio load/CA and the ratio 
Ship/Load to refine the estimation of the distribution for the ratio load/CA, 
which in turn improves the uncertainty of the forecast for the final CA. 
Figures 8 and 9 illustrate the situation of using this relationship. Figure 8 
shows the situation of forecasting for 2Q2003 CA for a particular computer 
product. On the graph, the dot labeled 3Q02 shows the actual for the ratio of 
load to CA. If this height of this dot (which is the ratio of load to CA) in this 
graph can be estimated precisely, then the forecast for the quarter CA for 
3Q2002 can also be estimated precisely (obviously, current load divided by 
this ratio gives the quarterly CA forecast). Without the fitted functional 
relationship from history, there is no other indication that would lead to the 
conclusion that 3Q02 load to CA ratio would be that low. Any model making 
use of history of this ratios would lead to a forecast of this ratio much higher 
than the actual (for example, the yellow line shows the weighted average of 
the ratios in history). This dramatically improves the forecast accuracy for this 
case. A similar situation is shown in Figure 9, where a comparison of forecast 
to actual for CA is made between BIA and a non-BIA method commonly 
used. 

Function block 501 in Figure 5 is the same as function blocks 31 and 
41 in Figures 3 and 4 for CA L and CA 5 , respectively. It creates the ratios of 
load to CA history collection. Function block 502 computes the corresponding 
ship-to-load ratio. Function blocks 503 to 507 compute the necessary 
quantities for determining the least square error estimate of the coefficients in 
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the equation = b — . Function blocks 508 and 509 use these 
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quantities to determine the parameters a, b. The derivations follow from 



setting up a series of equations like = b 
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for the zth history 



quarter and jth week. Take the log of both sides and form the sum of square 
5 error of both sides. Taking derivatives of this sum of square error and equating 

it to zero gives the condition for the parameters a, b to minimize the sum of 
square error. The results are closed form minimum least square solution as 
shown in function blocks 508 and 509. Function block 510 uses the estimated 
coefficients to determine the model fit error e, one for each historical data 

10 point. In function block 5 1 1 in Figure 5B, the sample sigma for such error is 
computed, which is also the sigma for the final forecast error. In function 
block 512, the computed least square solution of the model parameters a, b is 
used to determine the estimated load to CA ratio, which in turn is used in 
function block 513 to determine the CA forecast using current week quarter to 

15 date load L } by dividing the load with the forecasted load to CA ratio. This 

output 514 is the forecast of CA using Load/Ship (LS). 

The forecast generation from Load and CRAD performed in function 
block 212 of Figure 2 is shown in Figures 6A and 6B. This function in the 
invention exploits the relationship between the ratio load/CA and the 

20 statistical property of the CRAD for improving the forecast made by Load or 
Ship alone. The idea was triggered by the observation that when the orders are 
piling up but they are weak, the customer requested ship date (CRAD or RSD) 
generally is moved toward the end of the quarter just as a place holder. One 
reason this happens is that sales person often has a quota to make during the 
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quarter (before the quarter ends). If the demand is weak and they have 
difficulty making the quota, some of them would ask the customer to place an 
order just as a place holder. But because the customer is not genuinely 
interested in buying, the order generally cannot be scheduled to be shipped 
5 soon. As a result, the CRAD or RSD date is set more toward the end of the 

quarter. Without using such information, an artificially built up high order 
would naturally result in a high forecast. The idea here is to use the degree 
with which the CRAD histogram shifts toward the end of the quarter as an 
indicator to determine how much the normal forecast should be scaled back 

10 due to artificially built up order. One way to do this is to explore the 

relationship between the adjustment or offset for the historical quarters where 
the actual data are already known, so that the amount of adjustment needed to 
bring a forecast to actual is known, and the degree of shift of the CRAD's 
histogram. Figure 10 does exactly that. It is seen that the relationship is indeed 

15 present, with R-Square fit of 95.5%, and is log linear (if the SNR is the ratio 

of mean divided by the sigma of the CRAD histogram) or a linear function (if 
the SNR is the log of the ratio of mean divided by the sigma of the CRAD 
histogram), the adjustment needed can be predicted such that the forecast error 
is drastically reduced compared with the case where such adjustment is not 

20 made. Figure 1 1 shows exactly that. For a particular product family, it is 

shown that three quarters of history (represented by three dots, not counting 
the third one starting from the left) were used to estimate the adjustment 
model (shown by green curve), and this model is used to predict a new 
quarter's demand. The third dot counting from the left is the actual adjustment 

25 needed for the new quarter whose demand is to be forecasted. From Figure 1 1 , 

it can be seen that with the fitted model, the forecast error can be reduced by 
more than 35%, a huge amount. Figure 12 is another way to look at the system 
function of the CRAD method. 
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Now the method and the procedure shown in Figures 6A and 6B for 
this CRAD methodology will be explained. In function block 601 in Figure 
6A, the mean and sigma of the histogram are computed of all the orders with 
CRAD dates as the horizontal axis. This is done for every quarter i in history 
5 and every weeky. In function block 602, the SNR is computed. Note that the 

log is taken. When the log is taken, the regression model in Figure 10 becomes 
a linear rather than a log linear function. In function block 603, the adjustment 
needed to bring a normal forecast based on load to the actual for each of the 
historical quarter i and week j is computed. In function blocks 604 to 610 
10 (continuing to Figure 6B), the necessary quantities for determining the least 

square error estimate of the coefficients in the equation e .. = b +a . CA7J? are 

computed, where E, 7 is the adjustment on the normal forecast needed to take 
into consideration the weakness of the order as shown in function block 603. 
In function blocks 61 1 and 612 in Figure 6B, these quantities are used to 
15 determine the parameters a p b r The derivations follow from setting up a series 

of equations like e y . = bj+a jSNR _ for the zth history quarter and yth week. Form 

the sum of square error of both sides algebraically. Taking derivatives of this 
sum of square error and equating it to zero gives the condition for a p b } to 
satisfy to minimize the sum of square error. The results are the closed form 

20 minimum least square solution determined in function blocks 61 1 and 612. 

Function block 613 uses the estimated coefficients to determine the model fit 
error for the adjustment model 8, one for each historical data point. In function 
block 614, the sample sigma for such error is computed, which is also the 
sigma for the final forecast error. In function block 615, the SNR is computed 

25 for the current quarter using the current mean and sigma for the current CRAD 

distribution. Function block 616 uses the computed least square solution of the 
model parameters a p b } to determine the estimated adjustment needed, and add 
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it to the \i LCAj , which is the forecast made by using only the load and CA 
information in Figure 3. This completes the forecast of CA using Load and 
CRAD information. 

The adaptive optimization 218 in Figure 2 is shown in Figure 7. The 
function in Figure 7 is the last function of the process. Before the start of the 
function, seven forecasting outputs are directed to it. They are CA^, CA 5 , 
CA^, CA L CRAD , CA SCRAD , CA^c,^, and CA hisr The goal is create a final 
optimum forecast. There are several keys to this function. One is that it picks 
the best forecast model specific to the geography and product group to which 
the forecast is to be applied. This is very crucial to the success of the function. 
The second key is that it eliminates candidates depending on known properties 
according to how long the historical data are available and whether the time of 
forecast is early in the quarter, late in the quarter or not. 

A determination is made in decision block 71 as to whether a new 
quarter has just arrived and the actual for the old quarter just became 
available. If so, then decision block 71 will direct the process to go to function 
block 72 to update the forecast error performance metric e CAijk that is 
maintained for each geographic region j, each product group k and each week 
i. Decision block 73 will bypass the LS forecasting method in function block 
74 if the historical length is shorter than three or if the current week is still 
early in the quarter or very late in the quarter. Because the LS model fits a 
power regression with two parameters, it is essential to have at least three 
points of history so as not to overfit. Furthermore, when it is very early in the 
quarter, the ship is too small to make the LS work effectively. Similarly, 
decision block 75 will direct the system to bypass any forecast made in 
function block 76 with only ship if the week number is less than two, because 
the ship usually starts building up much later than load and is more prone to 
error. In function block 77, any method from the candidate lists based on any 
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information not available to the model is eliminated, based on human 
judgement. In function block 78, a search is made among the remaining 
candidates (for each geographic region, product grouping and for the current 
week) for the one that has the smallest mean average percent error based on 
weighted historical performances. The candidate that is picked will be the one 
chosen as the final forecast in output 79. 

Figure 12 summarizes the invention. Forecast generation 121 is first 
performed using CA hist 122, Load history 123 and current Load 124. Then, 
using the log mean to sigma ratio of the CRAD distribution 125, the CRAD 
history 126 and the current CRAD 127, the forecasts are refined at 128 to 
arrive at more accurate forecasts. Using adaptive optimization, a final optimal 
forecast is selected at 129. 

Best Indicator Adaptive (BIA) method is significant both in terms of 
theoretical foundations and practical impact and implications. The common 
theme and unifying theory of the power of quotient, and the methods of 
making use of order composition and sales opportunities pipeline progression 
as well as the methodology and theoretical analysis of the CA Quarter History 
indicator, and the adaptive optimization framework, are all key contributors. 

While the invention has been described in terms of a single preferred 
embodiment, those skilled in the art will recognize that the invention can be 
practiced with modification within the spirit and scope of the appended 
claims. 



